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Abstract 

We present a phenomenological study of the single-transverse spin asymmetry in azimuthal 
correlations of two jets produced nearly "back-to-back" in pp collisions at RHIC. We properly take 

> . 

, into account the initial- and final-state interactions of partons that can generate this asymmetry in 

' QCD hard-scattering. Using distribution functions fitted to the existing single-spin data, we make 

predictions for various weighted single-spin asymmetries in dijet correlations that are now readily 

In- 
testable at RHIC. 
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I. INTRODUCTION 



Single-transverse spin asymmetries (SSAs) play a fundamental role for our understanding 
of QCD in high-energy hadronic scattering. They may be obtained for reactions in, for 
example, lepton-proton or proton-proton scattering with one transversely polarized initial 
proton, by dividing the difference of the cross sections for the two settings of the transverse 



polarization by t 
asymmetries 



reir sum. There have been extensive experimental investigations of such 



. These have initiated much theoretical progress, in particular 
within the last few years. 

A particular focus has been on a class of single-spin observables that are characterized 
by a large momentum scale Q (for example, the virtuality of the photon in deeply-inelastic 
scattering (DIS)) and by a much smaller, but also measured, transverse momentum q±. In 
such a "two-scale" situation, single-spin asymmetries may arise at leading power, that is, not 
suppressed by an inverse power of Q. For some of these cases, factorization theorems have 



been established 



0,y, 



10] that allow to write the spin-dependent cross sections in terms of 



parton distribution functions and/or fragmentation functions, perturbative hard-scattering 
functions, and so-called soft factors. A crucial feature is that the distribution functions 
and the soft factor in this factorization are not integrated over the transverse momenta of 
partons, because these in fact generate the observed transverse momentum q±. Among other 
things, the observables may therefore provide valuable insights into the dependence of parton 
distributions in nucleons on transverse momentum. This becomes particularly interesting 
when the nucleon is transversely polarized, because there may be correlations between the 
nucleon spin vector, its momentum, and the parton's transverse momentum. One particular 
correlation, known as the "Sivers effect" and described by so-called "Sivers functions" 11( , is 
now widely believed to be involved in a variety of observed hadronic single-spin phenomena. 

Closer theoretical studies have revealed that the Sivers effect plays an important role 
in QCD, beyond giving rise to phenomenological functions to be used in the description of 
single-spin asymmetries. A particularly interesting feature is that the Sivers effect is not 
universal in the usual sense, that is, it is not represented by universal probability functions 
convoluted with partonic hard-scattering cross sections. This might at first sight appear 
to make the study of these functions less interesting. However, the non-universality has in 
fact a clear physical origin, and its closer investigation has turned out to be an extremely 
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important and productive development in QCD. In a nutshell, in order not to be forced 
to vanish because of the time-reversal symmetry of QCD, single-spin asymmetries require 
the presence of a strong-interaction phase. For the Sivers functions this phase originates 



from the "gauge links" in their definition [12|, 



14j, which are path-ordered exponentials 



of the gluon field that make the functions gauge-invariant. In DIS, the gauge link may be 
viewed as a rescattering of the parton in the color field of the nucleon remnant. That such a 
final-state rescattering may generate the phase required for a nonzero SSA in semi-inclusive 
hadron production in DIS (SIDIS) was first discovered within a model calculation [is] ]. 

Depending on the hard-scattering process, the "rescattering" will however manifest itself 
in different ways. For example, for Drell-Yan lepton pair production in hadronic scattering, 
initial- state, rather than final-state, interactions are relevant. As a result, the phase pro- 
vided by the g aug e links is opposite, and the Sivers functions for the Drell-Yan process have 
opposite sign [l2,[l^, li, 15]. In more general terms, the nontrivial "universality" property 



of the Sivers functions is the direct consequence of gauge interactions in Quantum Chromo- 



13 



dynamics I 12L ... , . ... 



14j, and of the QCD factorization theorems applying to the relevant hard 
processes y, LUJJ. It is a remarkable and fundamental QCD prediction that really tests 
all concepts we know of for analyzing hard-scattering reactions in strong interactions, and 
it awaits experimental verification. 

While measurements of SSAs in SIDIS are now maturing and have established the pres- 



ence of Sivers-type contributions [6|, it will still be a while until precise single-spin mea- 
surements in the relatively rare Drell-Yan process will become feasible at RHIC or 



elsewhere [TZj]. However, there are of course other hard-scattering reactions in hadronic 
collisions for which single-transverse spin asymmetries may be defined, and that may po- 
tentially be used in lieu of the Drell-Yan process for testing the nontrivial "universality" 
properties of the Sivers functions. In [181 ]. it was proposed to use the SSA in azimuthal cor- 
relations of two jets produced in pp collisions at RHIC to learn about the Sivers functions. 
To a first approximation, such dijets are produced by 2 — > 2 partonic QCD hard-scattering. 
With collinear kinematics, the jets are exactly "back-to-back" in the plane perpendicular 
to the initial beam directions and thus separated by 180° in azimuthal angle in this plane. 
Partonic transverse momenta will generate deviations from this topology, because they will 
lead on average to a non- vanishing net transverse momentum qj_ of the jet pair, much smaller 
than each of the jet transverse momenta P± individually. The observable is therefore of the 
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"two-scale" type described above, and as was shown in [18j, if one proton is polarized, a 
single-spin asymmetry may be defined that is in principle sensitive to the Sivers functions. 
As a caveat, factorization of the spin- dependent cross section for this observable in terms of 
transverse-momentum-dependent (TMD) functions still remains to be established. 

Unlike the relatively simple cases of SIDIS and the Drell-Yan process, where either final- 
state or initial-state interactions contribute to the Sivers asymmetry, both are present for 



dijet production 



19] . This complicates the analysis of the process-dependence of the Sivers 



functions considerably, but at the same time it also makes it much more interesting from 
a theoretical point of view, because the interactions in this case are "truly QCD", that 
is, they involve the detailed gauge structure of the theory, including for example its non- 



abelian nature. At the time of 



18| . the pro cess- dependence had not yet been worked out 



for the case of the SSA in dijet production, so that phenomenological predictions had to 
remain qualitative, at best. Over the last year, however, there has been extensive work on 
deriving and clarifying the structure of the gauge links for this and related processes in pp 
collisions l2ll . |22j . Indeed, the resulting structure is far more complicated than that in 
SIDIS or the Drell-Yan process. However, it turns out that if one takes a certain weighted 
integral ("moment") of the asymmetry, remarkable simplifications occur. This moment is 
defined by integrating the spin- dependent cross section with a factor sin 5, where S is the 
azimuthal imbalance between the two jets (5 = if the jets are exactly back-to-back in 
azimuth). For each of the various contributing 2^2 partonic channels, the gauge link 
then essentially collapses into a set of simple color factors that multiply contributions from 
color-gauge invariant subamplitudes to the given partonic process. One may, in fact, for 
convenience choose to absorb these factors into the hard-scattering functions, and define the 
Sivers functions as the functions measured in the SSA in SIDIS. In this way, the net effect of 
the gauge links on the sin <5-moment of the spin-dependent cross section is to generate new 
partonic hard-scattering functions that are different from the usual spin-independent ones, 
but that are actually similarly simple in structure. 

At the same time, taking the moment of the factorized spin-dependent cross section 
leads to a new expression that involves only a certain moment of the Sivers functions in 
partonic transverse momentum k±, rather than the fully transverse-momentum dependent 
functions themselves. As was shown in [l^j], these /cj_-moments of the Sivers functions are 



directly related to twist-three quark-gluon correlation functions first introduced in [231,124 to 
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describe the SSA for single-inclusive hadron production in hadronic scattering. By now, quite 
some knowledge about these correlation functions has been gathered from phenomenological 
studies [25j of the corresponding data. 

The upshot of all this is that it has now become possible for the first time to make 
predictions for the sin 5-moment of the single-transverse spin asymmetry in dijet production 
at RHIC that correctly take into account the process-dependence of the Sivers functions 
and incorporate phenomenological information on some properties of the functions that is 
available from other measurements. This is the goal of this note. 

II. SPIN-DEPENDENT CROSS SECTION AND sin 5-MOMENT 

We study azimuthal correlations of two jets produced nearly "back-to-back" in a hadronic 
collision. More specifically, we are interested in situations in which the sum of the two jet 
transverse momenta, q± = P\±+P2± (or a component or projection thereof), is measured and 
small, while both individual jet transverse momenta are large and similar. We will therefore 
approximate Pu_ = —P2± = P± wherever possible. For the lengths of these momentum 
vectors we will simply write P± = \P±\ and q± = \q±\. The differential cross section for the 
process with one transversely polarized hadron contains terms of the form 

+ e z ■ [S x x q ± ) - — - — — - , (1) 



drji dr/2 dP 2 d<pi d 2 q± drji dr\2 dP]_ d<pi d 2 qj_ 2 V / drji dr\2 dP]_ d(f>i d 2 q 

where e z is a unit vector in the direction of the polarized proton, S± is the transverse spin 
vector of the polarized proton, and rji and r^ 2 are the pseudo-rapidities of the two jets. 
The first term in Eq. (CQ) represents the unpolarized (or spin-averaged) cross section, while 
the second term is the single-transverse-spin dependent one. We note that the angular 
dependence of the spin-dependent term is |Sj_| sin(0fc — </>$), where <ph = (4>i + 4>2)/2 is the 
so-called bisector-angle of the two jets, which (approximately) corresponds to the direction 
of q±. For this reason one may also choose to integrate the six- fold differential cross section 
in Eq. ([T]) over 1; keeping q± fixed. 

As a generalization of the SIDIS and Drell-Yan cases |9j, we can write down a factorization 
formula for the differential cross section at small imbalance (q± ^ 0) of the jets, in terms 
of TMD parton distributions, soft factors, and hard-scattering functions 26|. We remind 



the reader that such a factorization still remains to be proven. In this paper, we will not 
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discuss the details of factorization issues related to the dijet correlations. As we mentioned 
above, significant simplifications occur when the imbalance of the two jets is integrated out 
by taking certain moments. For example, integrating the spin-independent differential cross 
section over all q±, its expression reverts to the standard collinear factorization formula for 
dijet production, 

= f , 2 27r dPauu 27T d A a vv 

UU J ± drji dr/2 dPj_ d(j>x d 2 q± dr\\ dr\2 dP]_ d(j>i 

= ^2x a f a (x a )x b f b (x b )H^ cd (s,iu) , (2) 

ab 

where the /, denote the usual collinear (light-cone) parton distribution functions for par- 
ton type % = q, q,g. We have assumed here that these are recovered by integration of the 
corresponding TMD distribution functions over all partonic transverse momentum, and we 
disregard the renormalization properties of the operators defining these distributions and 
the scale dependence introduced by renormalization. However, all these effects can be sys- 
tematically included accordingly. 

The factors H^ b ^ cd in Eq. (J2]) are the customary hard-scattering cross sections do ah ~' cd j dt 
for the partonic processes ab — > cd (for a compilation, see, for example, Ref. [22|). They are 
functions of the partonic Mandelstam variables, s = (p a + Pb) 2 , t = (p a —Pc) 2 , u = (p a —Pd) 2 , 
with obvious notation of the parton momenta. In terms of the hadronic center-of-mass 
energy a/s and the jet transverse momenta and pseudo-rapidities, one has s = x a x b s, t = 
—P 2 (e 712 ' 711 + 1), u = —P 2 (e m ~ V2 + 1), where the partonic momentum fractions are fixed 
by x a = £k (e" 1 + e" 2 ), x b = ^ (e"" 1 + e"^). 

Next we turn to the single-transverse-spin dependent differential cross section daru i n 
Eq. (TjQ), which can be further simplified by taking a moment in 

d . c e z -(P ± x q ± ) 

P± sin 5 = , (3) 

-fi 

where 5 = it — (02 — 4>i) measures how far the two jets are away from the back-to-back 
configuration. Within our approximations, the weight P± sin 5 corresponds to a weight in 



q±. One finds [211 ] : 



2 Px sin S \ i ri i f j2 -* 9-L 27r d G axu 



Mp I ' TU ^ ± J ^ X Mp drjidrj 2 dP'j_d(l)id' 2 q 1 



S±\^2x a — gT Fa (x a ) x b f b (x b ) HJ™ s cd (s, i, u) , (4) 



ab 



where Mp is the proton mass, g is the strong coupling constant and the Tp a (x) are the 
Qiu-Sterman matrix elements or quark-gluon correlation functions 

d£~dr]~ 



T Fa {x) 



An 



241 ]. defined as 

''^^S^(PSmoh + F a ^r)r{v-)\P3) , (5) 



with the quark fields if) a (for flavor a) and the gluon field strength tensor F a + . The Tp a (x) 
matrix elements enter because they are related to the fcj_-moment of the (SIDIS) Sivers 
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function for quark flavor a [141 ]. that is, gTp a {x) = -2M/ 1T(1 (i). We note that there could 
also be contributions by purely gluonic u ggg" correlation functions. These will be ignored 
for now, so that the label a in Eq. (jlj) runs only over quarks and anti-quarks. Furthermore, 
there is actually a second contribution to the single-transverse-spin dependent cross section, 
which involves the scattering of transversely polarized quarks from both the polarized proton 
(transversity distribution. 8f or hi) and from the unpolarizcd proton (Bocr-Muldcrs function 
hi 23]). The latter functions are also effects of initial- and final-state interactions and appear 
in a matrix element for unpolarized protons similar to Eq. ([5]). Like the Sivers- type u ggg" 
correlation functions, we will also ignore the contributions associated with the Boer-Mulders 
functions in the present study, even though a future more detailed analysis of forthcoming 
experimental data may well require to take all of these into account. 

The relevant hard-scattering functions H^^ d have been calculated in 2l|, l22| , where they 
were termed "gluonic-pole cross sections" due to their association with the Qiu-Sterman or 



gluonic pole matrix elements 



241 ] . We have also independently reproduced [z 



the # a6 _> cd 



within a model calculation along the lines of Ref. [15(. For convenience, we list the ones that 
will be relevant for our numerical calculations presented below: 

a 2 n N 2 - 5 2(s 2 + u 2 ) 



rrsivers 
qq'^qq 



rrsivers 
n qq'^qq 



,(s,t,u) 
f(s,i,u) 



a 2 ir 



AN 2 
N 2 „ 



i 2 

3\ 2(s 2 + u 2 ) 



Hqq™q'q-'\Sj ^, It) 



Hqq-^qq(.Si ^? ^) 



AN 2 



a 2 s iiN 2 + l2(t 2 



t 2 



u 



s 2 AN 2 
a?vr f N 2 - 5 



a 2 ir 



AN 2 



2(s 2 + u 2 



N 2 -3 2(s 2 



u 



AN 2 



+ 



+ 



2(s 2 + t 2 



N 2 + 3 4s 2 
AN? 



N 2 + 1 2{u 2 + i 2 ) N 2 c +lAu 2 



AN 2 



AN? St 
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N 2 



2(s 2 + u 2 



4(N 2 - 1) 



t 2 

2(s 2 + u 2 



u 



rrsivers , 



S, t, It) 



a 2 ir 



N c 2{t 2 



AN 2 (N 2 - 1) su 
1 2(t 2 + it 2 ) _ 
'2N C s 2 ^4 
2N 2 + 12(t 2 + it 2 



u 



2(s 2 + it 2 ) 



2(N, 



t 2 



t u 
- + - 
u t 



(6) 



4iV c 3 tix 

where N c = 3 is the number of colors. Hard-scattering functions corresponding to gluonic 
correlation functions were also calculated in [2^], but are not taken into account in our 
present study as we discussed above. Comparing the above functions with the usual spin- 

n 

averaged hard-scattering functions [22], one can see that they essentially differ in the color 
factors that multiply terms with similar kinematic structure. This is the net result of the 
combined effects from the initial and final state interactions in dijet production in hadronic 
reactions. 



III. PHENOMENOLOGICAL RESULTS FOR RHIC 



We now use the above formulas to obtain some predictions for the SSA in dijet-pro duct ion 
at RHIC. A SSA for our (2P± sin <5/Mp)-moment (i.e. gj_/Mp-moment) can be defined as 
(from now on we choose |Sj_| = 1) 

(2P ± sin 5 /M P ) TU _ Eab x a ^gT F a {x a )x b f h {x b )H^l d {s, i, u) 

( X )UU Y,ab X afa{x a )x b f h (x b )H2^ cd {s i tAL) 

The quark-gluon correlation functions T Fa have recently been fitted [25j to data |2j, |5| for 



the SSA in sing 
presented in 



e-inclusive hadron production in hadronic collisions. Two sets of T Fa were 
25]. For definiteness, we choose set I, which is a two-flavor fit with valence u 
and ci-quark T F distributions only. For these the following parameterizations were obtained 
in [25|: 



T, 



Fu v \ 



0.275a; 



0.508, 



1 — X 



,0.399 



u v {x 



Fd v {X) 



-0.365a; 



-0.108/ 



(l-a;) - 287 ^(a;) , (8) 



where the u v ,d v are the corresponding unpolarized valence quark distributions. For the 
latter, as for all other unpolarized parton distributions we need, we choose the CTEQ5L 
parameterizations 29| . We plot the resulting weighted asymmetry in Fig. [TJ The kinemat- 
ics are chosen to correspond to current measurements at STAR: both jets have transverse 
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momenta 5 GeV < P± < lOGeV and pseudo-rapidities —1 < r/i < 2. The asymmetry is 
plotted as a function of the sum of the two jet pseudo-rapidities. We have chosen the hard 
scale /i = P± in the strong coupling constant and the unpolarized parton distributions. 

For comparison, we also show in Fig. [T] the asymmetries that one would have expected 
if the Sivers functions relevant for dijet production were identical to the ones in SIDIS or 
the Drell-Yan process. In the framework of our present calculation, the corresponding par- 
tonic hard-scattering functions would in this case be identical to the spin-averaged functions 
Hab-^cd' or t° their negatives. The dotted-dashed lines in the right panel of Fi g. [1] represent 
these cases. They essentially correspond to the earlier predictions of [3] and 261] that were 
made when the process-dependence of the Sivers functions had not yet been worked out for 
the dijet case. One can see that when the correct pro cess- dependence is incorporated, the 
asymmetry overall becomes smaller by about a factor two, which can be traced back to the 
color factors for the dominant subprocess qg — > qg. The sign is identical to the case when 
the Sivers functions for dijet production are assumed to be "DIS-like", implying that in a 
sense final-state interactions dominate over initial-state ones. 

In principle, one might verify experimentally the process-dependence of the Sivers func- 
tions by discriminating between the various curves in Fig. [1] and confirming the QCD- 
predicted result shown by the solid line. In practice, this may require good knowledge of 
the Tp distributions, and an understanding of issues like scale evolution and higher-order 
corrections. A closer analysis reveals that the asymmetry in Fig. [T] is the result of rather 
significant cancellations between contributions of opposite signs by u and d quarks. To show 
this, we display their individual contributions separately in the left panel of Fig. [TJ 

We finally also briefly discuss a related type of SSA in dijet production at RHIC. In [261 ]. 
a differently weighted SSA was considered, defined in the following way: 

where Sgn(x) is the sign function. This asymmetry has the property that the weight only 
depends on the azimuthal separation of the two jets, but not on their transverse momenta. 
This may be an advantage for experimental measurements when the jet energy scale is not 
precisely known. When applied to the TMD factorized expression for the spin-dependent 
cross section, the moment defined in Eq. leads to an expression different from One 
finds in fact that in general the transverse-momentum dependences of the various functions 
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0.1^ DijetSSA: 5GeV<P i <10GeV, -1<n i2 <2 



(2P i sin8/M p -weighted) 



0.1^ DijetSSA: 5GeV<P i <10GeV, -1<H 12 <2 



(2P i sinS/M p -weighted) 



0.05 



0.05 



-0.05 



-0.05 



y\ 4 



yii 2 4 
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FIG. 1: The weighted single-transverse-spin asymmetry for dijet- correlations in polarized proton- 
proton scattering at RHIC, as a function of the sum of the jet pseudo-rapidities. The jet transverse 
momenta are in the range 5 — 10 GeV. The solid line is our main result, based on Eqs. fiji"^. In 
the left panel we also plot the individual contributions due to the u-quark (dashed line) and d-quark 
(dotted line) Qiu-Sterman matrix elements. For comparison, in the right panel we also show the 
results that would be obtained if the Sivers functions contributing to dijet production were the same 
as those in SIDIS (upper curve) or the Drell-Yan process (lower curve). 



do not completely decouple anymore, so that the final result can in general no longer be 
expressed in terms of only functions of light-cone momentum fractions. If one assumes for 



simplicity 



261 ] , however, that the only relevant dependence on transverse momentum resides 



in the Sivers functions, the resulting expression again resembles a collinearly-factorized one: 
(Sgn(5)) TC/ Z ab ^Ta ] (x a )x b f b (x b )H^T cd (s, i, u) 



(1> 



uu 



Eafe Xafa{Xa)X b f b {x b )H2^ cd {s, t, u) 



(10) 



where the q^l^ are the so-called "1/2-moments" of the Sivers functions for SIDIS and are 
defined as 



d 2 k\ 



M 



(11) 



In the above equations, we have used the notation "q^a" for the Sivers function for quark 
flavor a; its definition is identical to that in [ijj]: qx a = frr a - m 26 ] , the hard-scattering 
functions were chosen to be the same as the spin-averaged ones, with opposite sign. In light 
of the above discussions, one now would like to update the predictions for the asymmetry 



10 



0.1^ DijetSSA: 5GeV<P i <10GeV, -1<n i2 <2 



0.1^ DijetSSA: 5GeV<P i <10GeV, -1<n i2 <2 
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0.05 
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-0.05 



y\ 4 



■« u « T|+T| 

'2 '1 '2 

FIG. 2: Same as Fig. 1, but for the SSA defined in Eq. \10\) . 



in (fTUj) . We will use the H a ^^ d given in Eqs. ([HD- We note that it remains to be established 
that the same Hf^* d do contribute in this case. This is not a priori clear, because the 
general gauge link structure is very complex in the general TMD case, and it has so far only 
been demonstrated that the H^^d apply when the sin 5 moment is taken. For now, we just 



conjecture that use of the H a ™^ * d of [2ll . |22J is justified in this case; a closer discussion of 



this issue is left for future work 



28] 



The 1/2- moments of the quark Sivers functions were determined in 26] by a fit to the 
HERMES data |6j for the Sivers-type single-spin asymmetry: 



u 



(1/2), 



— 0.75x(l — x)u(x), d T 



(1/2), 



X) 



2.76x(l - x)d(x) 



(12) 



These results correspond to set II presented in 26j. In Fig. [2J we show predictions for the 
asymmetry for the dijet-correlation defined in (flOl) at RHIC, based on this parameterization. 
We find that the asymmetry shares many features with the one shown for the (2Pj_ sin 5/Mp)- 
moment in Fig. [TJ We note that first preliminary experimental data for this asymmetry have 
now been presented by STAR 3], which are so far consistent with a vanishing asymmetry. 



IV. CONCLUSIONS 



In summary, we have studied single-transverse spin asymmetries in dijet-correlations 
at RHIC, making use of the recently derived partonic hard-scattering cross sections that 
properly incorporate the initial- and final-state interactions, and of distribution functions 

11 



fitted to existing data for single-spin asymmetries. We have found that the initial- and final- 
state interactions tend to decrease the magnitude of the asymmetry at RHIC with respect 
to earlier estimates that assumed the Sivers functions for this observable to be identical to 
the Sivers functions in SIDIS or the Drell-Yan process. Overall, the resulting asymmetries 
turn out to be more dominated by the final-state interactions, and hence more "DIS-like" . 

With experimental data on dijet single-spin asymmetries now forthcoming, it will be 
interesting to perform detailed comparisons with the theoretical expectations. Other ob- 
servables, such as the SSAs in dihadron production pp — > hih 2 X or in photon-plus-jet 
production pp — > 7jetX, will also be extremely interesting. It will be important to further 
develop the theoretical framework for all these observables, by addressing issues like TMD 
factorization, higher orders, soft factors, and Sudakov suppression, in particular. 
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